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ABSTRACT 


A computer data reduction system has been reprogrammed to 
sample individual scan recordings of the line-spread function 
of laser light propagating through the turbulent marine boud- 
ary layer. 

The distribution of spatial frequencies of the Modulation 
Transfer Functions from random single scan samples from 
experimental data sets has been investigated. A preliminary 
analysis indicates that the zero spatial frequency component 


of the MIF is log-normal distributed. 


oe 
a 


of sere 
a 


sy ato’ dash rate he 


aT ; : 
—i we > ~~ ee p 
Ho yettegaqote tdgtl ae 
i 


ui. 


Tit. 


TABLE OF CONTENTS 


ENTRODUCTION: xian td atcake th Ge Sere oom ed eds owes 9 
ANALYSIS“ O8 DAWA: 126 eh tnao each bESGee Rieke eka Rae et 11 
A. (FUNCTIONAL DESCREPTION ..cui2cnekanwa ne deeea¥er 11 
De “Deleas gcse Sees pauy Seon we eb es i eaeAar mt 
2, Data Reguction 200440242408 ennces Red wie es it 
S.. Statistical Analysis. .coccusudtestedakeek Ss 12 
B. PROGRAM THEORETICAL CONSIDERATIONS ............ 18 
i, (NOM Cube’ oe & oe ee oes be ee So teens 18 
2. “GOSdneSs. Of TIC «22 sbyed soe cuenta eeaceesds Lg 
C. COMPUTER PROGRAMS ........ cc eee ee eee eee te eee 20 
1, Nie 80 Data Reduction Proeram .....6se655%5 20 
2, Data Analysis Programs: «icc¢sdi4eeegutowean 20 
5. Program Ploweneres: «0134 das u8¢eaSeeeseeeees Z2 
RESULES © ct ogi dee oid dea aa ante, @ wae te aay Biv or ao ee aad ee eed we 30 
A CONSISTENCY OF DISTRIBUTION ................... 30 
B FACTORS AFFECTING DISTRIBUTIONS ............... 5) 
Let Sern t al rOn: ok 26 ooh owe See ce tae ou keee a's 35 
2, ‘Sipmal @o Neise- FAC O- 4 5.neanna tees See aoe oo 
C. RESULTS AND DISCUSSION ...................220005 a0 
L. Log Normal Distribution of 
O Spatial Brequeney «to vawe dar deed ap be kets Sib) 
2. Normal Distribution of 
Higher Prequenties: oi.) «s dessa de eeades eeas 36 
3. - DiStribution versus Ca Value. .<es.e ene chases 36 
D. RECOMMENDATIONS FOR FURTHER WORK .............. 37 


STHeT HOD Se -GURAT 


ee | 
ro ; Tad QO. 


“abe ; VOMVeriozas J 
=), 
he 
ea. , 
“4 nat) obbed” «304 
‘ey par : ind setsase 
nyt od 
r sq “th 
a! 
; « “a 
ag 
SLA yere Fo Sarat 
ITUGIMTELE OUMTISTLA BSCR oo 
nolvai Jisalge 
S7n8 pelo 7c Jemgie oy 
‘ : 


HOTSEISSIG ame Whe 


; niatswcbsteh somabe Peek 
Page oe 1 is 


Or ee, & 


APPENDIX A: Operating Instructions for Data Analysis ...113 
COMPUTER PROGRAMS LISTINGS 


l. Data Analysis Program Listing ............e00000. 117 
2. Data Summary Program Listing .................... 120 
a. Mie 80 Program. Listine :2..0sSeccckaue cud see kwe as 121 
LUST OF REFERENCES cicsa 029 6640 e 4 SRR SRE RERHE SOAS Se HRS RO 125 
INITIAL DISTRIBUTION LIST wns cccas cee nee ie ete eee eee eet LZ6 


ee a ‘ : 
: i i af 7 ‘ ; may 


\ ‘ rare ‘ ‘ be 4 
Cif... OpewesmA §3a0 tor PAetPousTe- 


i 
oO 


Oo dG N WD wo fF Ww NYP FHF 


LIST OF DRAWINGS 


Sample Date SUmMary v4.06 sees etka di ems) Sess eho a we aS 14 
Sample MIF + 1 Standard Deviation PLOt ..40 060000 sees 15 
Sample Relative Distribution of MUP .2<viceseeteua ves 16 
Typical, Line: Spread Funetion. v.24 g105 644s teu dt eeweeades 17 
As 5B PrOeram OreaniZation: 140+ id«vevswueeddecicase ees 29 
A= D Data Analysis Program: 2.6 168 2634946645 tees ES 24 
A = 5 Data Summary Prostam uss e<ss ios PbS We wks Pee se 28 
A= G Data Analysis Program Listing «<siisavedioegcies 117 
Data. Summary Prosran Listing «. 40s o eens eee ered esagas 120 
A= De Nie 80 Program Listine dee eek st sd etendban dunes £21 


intel Ae Ahk - As 


~ rag -_ ar 
atta ao ter.” 


soltc. inwiyye “boabese? 0b ae 
| ey: 
= y shale Ss; 20 eS ty 
ie < boesae ous d 


Noting tao x0) vara 


ae 


s - 


sergo7? Stay leet al 
is hee aveye uc 7 


= 
ron at 


IV. 


LIST OF TABLES 


Data ANAlvsiS Sue ry: 64 ete nee we Pew oe ee eee Sad 31 


Grading Consistency of Data from 
Hae Dale SEU ¢ckene deine dk Leen stk semaey 6e40ed ee yawns 34 


Consistency of Individual Frequency 
Components from Samples of Same 
Data SE. 235645 28odeu £45554 Ss Se ese ede CRs ke GE AS eG 34 


Ratio of Distribution Graded 1 for 
Dit erent. (hn. Ranees cv ccee bakes H4bUen eee wee kawe vs 37 


at - oun waned ¢ dW Bev d & 1 OF = _esnene. 


= " xarsaper’ Losblitead My 
smae 20 caJgaue ac 
a tee t wamea hie @ 6 » ° rhe teen! 


3023 by hax0 ne Mutha be 
Se = 15 é é ‘ + *b ‘ ~as Sempaas co 


( 


I. INTRODUCTION 


Atmospheric turbulence adversely affects the propagation 
of a laser beam or the performance of an imaging system. The 
Modulation Transfer Function (MIF) is a measure of the degra- 
dation of resolution and thus the relative ability of the 
atmosphere to transmit spatial frequency information. Knowl- 
edge of the MIF of the atmosphere is necessary for prediction 
of the beam cross-sectional profile on target and spot wander 
for a high energy laser and its beam projection system. 

Members of the faculty of the Naval Postgraduate School 
(NPS), Department of Physics and Chemistry, have investigated 
the effects of atmospheric turbulence on the propagation of 
optical radiation in the marine boundary layer. 

The work started in 1973 and has continued through 1978. 
The research has been concentrated on developing a model for 
the MIF depending on wavelength, range and a properly path- 
position weighted value of the turbulence structure constant 
for optical index Cn. The work for the period 1 July 1974 - 
30 September 1977 is described in Ref. 1. Field experiments 
were performed by transmitting laser light of different wave- 
lengths along various seapaths. The received intensity was 
detected by a slit-scanning telescope system and recorded on 
tape. The Cn value was determined by computer calculations. 

The objective of this report is to describe the develop- 


ment of a computer program which samples individual scan 
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recordings of the line-spread function and determines a pre- 
liminary statistical distribution of the spatial frequency 


content of the Modulation Transfer Function. 
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II. ANALYSIS OF DATA 


A. FUNCTIONAL DESCRIPTION 
l. Data 

Experimental data are the intensities of the lLine- 
spread functions as recorded on a 4-channel magnetic tape. 
One track is used for verbal identification of date, time 
and pertinent experimental information. Two tracks are used 
for recording scans and the last track is used for recording 
the trigger signal. 

Verbal information has been cross-checked against 
experiment log book and MTF data output. Data have been 
recorded for one and two wavelengths simultaneously. Only 
data for .6328 microns have been analyzed in this report. 

The line-spread function was sampled by the scanning detector 
at a rate of 50Hz and recorded at a tape speed of 37.5 inches 
per second. 

2. Data Reduction 

The data reduction is done on a Nicolet Instrument 
Company, NIC-80 data processor. This is a data digitizer 
followed by a 12K 20-bit digital computer and a Tektronix CRT 
display unit. The NIC-80 computer is controlled by a Hewlett- 
Packard HP-9825 desk top type programmable calculator. 

The NIC-80 computer reads the line-spread function 
for one scan, performs a Fourier Transform (FT) analysis of 
the signal and transfers the FT spatial frequency components 


0 - 19 to the HP-9825. The MIF of the atmosphere is calculated 
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by dividing the FT of the line-spread function by the FT of 
the optics. The FT of the optics is stored by the NIC-80 and 
transfered to HP-9825. The experimental analog signal is 
digitized at 1024 points. The DC component is removed by 
manual setting of signal bias based upon observation of scan 
on CRT display. Any remaining DC component is compensated by 
subtracting the average intensity value of 32 points at each 
end of the digitized signal. 

When not performing FT and transfer operations, the 
NIC-80 remains in signal display mode with successive scans 
being displayed on CRT. 

The computation and transfer time for one scan has 
been observed to be 5 seconds, and will be the time between 
sampling of scans. During data reduction, the tape recorder 
was operated at 3.75 IPS, corresponding to 5 scans per second. 
This implies that about every 25th scan of the recorded data 
was sampled for analysis. The number of scans sampled was 
limited by the length of experimental data recorded, normally 
in the order of one minute, or 3000 scans. This allowed a 
sample size of 100 - 125 for data analysis. An attempt to 
run the recorder at .375 IPS speed was not successful due to 
problems with the trigger signal. 

3. Statistical Analysis 

The HP-9825 calculates the MIF at 20 different spatial 
frequencies from 0 to 38000 per radian in steps of 2000 per 
radian. These spatial frequency points are labeled 0 to 19 


and are tested for normal and log-normal distribution. Goodness 
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of fit is determined by the Chi Square test. A data summary 
containing mean, standard deviation and Chi Square values 
for normal and log-normal tests is outputted on the printer. 
Mean and + 1 standard deviation for frequency points 0 - 19 
is plotted, with the zero spatial frequency normalized to a 
value of unity and all other data adjusted correspondingly. 
A histogram showing relative distribution of MIF values for 
one specific spatial frequency is plotted as requested in 
program input. 

Examples of the above are given in Figures 1, 2 and 3 
respectively. The relationship of the MIF distribution to 
the Cn value has been investigated. The Cn values ranged in 


=O) 66 een 


the order of 8x10 Typical line-spread functions 
for two different Cn values are shown in Figure 4. Note that 
Cn squared values are given. 

The value of Cn used has been taken from existing MIF 
data outputs. If no output was available, an MIF analysis 
was performed. A detailed discussion of MIF and Cn is given 
in Ref. 1. The scan presentation is a function of time and 
detector position. The dwell time setting was obtained by 
analyzing calibration data recorded on tape. Dwell time is 
the sampling time per data point. Data have been obtained 
by recording the detector output with a diffraction grating 
of known dimensions placed in front of the telescope. This 
provides a known reference of angular spread. The dwell time 
is adjusted so that the angular spread of the first order 


diffraction as observed on the display is equal to the ref- 


erence value. 
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B. PROGRAM THEORETICAL CONSIDERATIONS 
The data were tested for normal and ln normal distribution. 
Ll. Normal Curve 


The normal or Gaussian curve is expressed as 
pe 
Y = constant exp arr ) 


where 


x is the value of MIF for normal test distribution 
or value of In MIF for log-normal distribution. 


x is the unbiased estimate of the Mean of the MTF 
values or Mean of In MTF values. 


s is the unbiased estimate of the standard deviation 
of the distribution. 


Y is the probability density of the distribution. 


Constant = i/s/¥27, where i is the width of a histogram 
bin to make area under curve equal to i. 


The Mean x is the arithmetic Mean given by: 
x = 


where N is the number of observations (samples). 
The standard deviation s is the positive root of the 


Variance s squared given by 
ee! 
(x, -%) “/ (N-1) 
a 


For computational purposes the following algebraically 
equivalent formula has been used: 
2 N N 


s- = ( 5 (x,)*-( 5 x,)7/N)/(N-1) 
i=l * i=1 + 


18 


The probability density curve has been computed and 

plotted over a range of Mean + 3 standard deviations. 
2. Goodness of Fit 

Test of goodness of fit to tested distribution has 
been done by histogram and Chi Square calculation. 

The width of the histogram bins is a program variable 
expressed as a ratio of the standard deviation. 1/3 has been 
used corresponding to a total of 18 bins. The number of 
samples with values within a bin is expressed as a fraction 
of 1, so that the area of the histogram is equal to the width 
of the bin interval. 

The Chi Square goodness of fit value has been calcu- 


lated from 


M (N£(n)-F(n)]7 


x > 
i<j, “OD 
where 
N = number of observations 
n = number of histogram bins 
F(n) = number of observations within a bin 


The term Nf(n) is the theoretical number of observations 


which should fall within a bin. The probability area has been 
calculated as the sum of a rectangle and a triangle approxi- 
Mating the area under the normal curve. The sampling 
distribution of the Chi Square statistics is tabulated as a 
function of degrees of freedom and percentile. The number of 
degrees of freedom is calculated from 


df = k-n-1 
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where 


= 
Il 


number of bins in histogram 

n = number of parameters estimated 
For a 18 bin histogram with mean and standard deviation 
estimated, the degrees of freedom is 15, with a 95 percentile 
value tabulated as 25. A calculated Chi Square value less 
than 25 will, at the 5% significance level, imply that the 


data fits the given distribution. 


C. COMPUTER PROGRAMS 
l. NIC-80 Data Reduction Program 
This program has been developed utilizing subroutines 
from the existing MIF analysis program. Pages 3 and 4 of the 
MIF program have not been changed. Main program is on page l 
and a data transfer subroutine has been inserted on page 2. 
Program flowchart is given in Figure 5A. Program listing is 
on page 122. 
2. Data Analysis Program 
Due to limited HP-9825 memory capacity, this program 
has been split into two parts stored on separate program 
files. During operation, programs are loaded and executed 
under program control. Program "A" is on file 8 and contains 
NIC-80 control and all analysis functions except data summary, 
which is on program ''B'' located on file 9. Programs are 
documented on flowcharts (Figures 6 and 7). Program listings 
are on pages 117-121. Execution of analysis program requires 


HP tape inserted with track 1 available for data storage. 
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Each data file must be minimum 160 bytes long and the number 
of files must be equal to or larger than the number of scans 
sampled. 

Execution time for a 100 scan analysis run is of the 
order of 90 minutes. Most of this time is required for tape 


movements and operations. 
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Figure 5A. Program Organization 
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Figure 5B. Program Organization 
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Date: 6 figure #, date, month, year 
Time: 6 figure # 
Cn: As argument and power 


Line 18 
Calculate 
Mean & For MIF spatial frequencies 0-19 
Standard Line 19-21 


Deviation 


Print , : 
MIF 0 & ae ee absolute and relative 

me Line 22 

Enter Histogram class interval is expressed 
Class as a ratio of standard deviation, e.g. 3 


Interval / Line 23 


Subroutine which plots 
' Mean plus-minus 1 standard 


' W 
ieee | deviation for spat. freq. 1-19 
Line 85-104 
NO 
YES To subroutine "LN" where natural 
ie logarithm of data is taken 
Line 76-84 


YNO, or already LN data on file 


Figure 6A. Data Analysis Program 
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Data 
Summary 
Program 
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Figure 7 


Enter Spatial frequencies 0-19 


Spat. Line 28 
Frequency 
Calculate X Scale equals mean plus/minus 
Plot 3 standard deviations 
Scales Line 30 


Relative distribution of MTF 
Line 32-42 


Calculates and plots normal curve 
based upon sample estimate of mean 
and standard deviation 

Line 43-46 


Figure 6B. Data Analysis Program 
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Figure 6C. Data Analysis Program 
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Figure 6D. Data Analysis Program 
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Figure 7A. Data Summary Program 
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III. RESULTS 


A total of 24 data analysis runs have been performed. 
Five of these have been on the same data sets in order to 
investigate consistency of distribution within the same data 
set. 

An arbitrary grading system for goodness of fit over a 
range of spatial frequencies has been established as follows: 


dig 


- All components in range within 95 percentile 


- One component in range outside 95 percentile 2. 


- Two or more components in range outside 95 percentile 


Data summary and MIF mean and standard deviation plots 
are given on pages 39-86, arranged in order of decreasing 
order of Cn. A distribution analysis summary is tabulated 
on pages 31-33. A set of histograms for all spatial frequen- 


cies for one data set is given on pages 87-112. 


A. CONSISTENCY OF DISTRIBUTION 
Data set number 12 has been run 5 times to verify consis- 
tency of distribution as only every 25th scan was sampled. 
Utilizing the grading system for spatial frequency ranges 
0, 1-6, 7-12 and 13-19, the following results as given in 


Table II and III, were found: 
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TABLE I 
DISTIRUBTION ANALYSIS SUMMARY 
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1.07| 031176 
Loe? 


O92177 


‘tio 
G25) 121776 
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Legend: 1. All values within Chi Square 95 percentile. 
2. One value above Chi Square 95 percentile. 
3. Two or more values above Chi Square 
95 percentile. 
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Legend: 1. All values within Chi Square 95 percentile. 
2. One value above Chi Square 95 percentile. 
3 Two or more values above Chi Square 
95 percentile. 
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TABLE II 


Grading Consistency of Data from Same Data Set 


Number of Runs with Grading 


Frequency Range Normal Log Normal 
Grade Lt. 2.2 LZ. 3 
0) 4 0 4 1 0 
1 - 6 1 0 4 >. *0.- Q 
p= 0 0 5 3. & 
13 - 19 GO @. 3 3 2° 0 


The number of samples which met the 95 percentile criterion 
for each spatial frequency in the 0 - 6 frequency range is 


listed in Table III. 


TABLE III 


Consistency of Individual Frequency Components 
from Samples of Same Data Set 


Number of Samples within 95% 


Frequency Normal Log Normal 
0 i 4 
i 
2 2 5) 
3 3 zs) 
4 3 5 
5 Z 5 
6 3 5 
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The results indicate that the distribution can be expected to 
give same Chi Square values for different sample series from 


same data set. 


B. FACTORS AFFECTING DISTRIBUTIONS 
1. Seintillation 
If the intensity of a laser beam which has propagated 
through a turbulent medium is observed at a point, the inten- 
sity appears to fluctuate or scintillate in a random way. 
It has been shown that the values of the logarithm 
of the amplitude follow a Gaussian distribution (Ref. 1, 
page 15). Spatial frequency component zero represents the 
area under the curve (signal). This area is proportional to 
the received intensity. It is thus expected that the distri- 
bution of the zero spatial frequency component should fit 
the log normal distribution. 
2. Signal to Noise Ratio 
The amplitude of the spatial frequencies decreases 
with increasing frequency number. For spatial frequency 15 
and above, the amplitude is of the order of 10% of the zero 
spatial frequency component. At higher frequencies noise 
will contribute more relative to the amplitude. Noise is 
considered normal distributed, and it is expected that the 


distribution at higher frequencies is normal. 


C. RESULTS AND DISCUSSION 
1. Log-Normal Distribution of 0 Spatial Frequency 
In 13 of 20 runs the spatial frequency zero fits the 
In normal Chi Square 95 percentile criterion. For runs which 
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do not meet the distribution criterion, the ratio of mean to 
standard deviation is in the order of 2. This may indicate 
saturation, and thus not follow log normal distribution (Ref. 
1, page 16). 
The results do indicate that the zero spatial freq- 
uency is log normal distributed. 
2. Normal Distribution of Higher Frequencies 
For data analysis the modulus, e.g., absolute value 
of the data has been used, and sign of the data has not been 
included. This will introduce a skewness of distribution. 
It is believed that this is the reason for lack of normal 
distribution for data at higher frequencies. 
3. Distribution versus Cn Value 
Data have been arbitrarily grouped into ranges with 
respect to Cn values as follows: 


Cn High: Cn > 6.107° 


8 8 


Cn Medium: Cn between 2.10 ~ and 6.10 


Cn Low: Cn < 2.1078 
The ratio of distributions graded 1 for the various 
Cn ranges is listed in Table IV. 


The results do not indicate that the distribution 


does change with Cn values. For all analysis results, it 


should be noted that the sample size is small, and the results 


should be interpreted as indications only. 
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TABLE IV 


Ratio of Distribution Graded 
1 for Different Cn Ranges 


Spatial 
Cn Frequency Normal Log Normal 
H 0 5/6 5/6 
I Loe * 2/6 3/6 
7 = 12 3/6 2/6 
H 13 = 19 0/6 0/6 
M 0 af? 4/7 
: 1 - 6 Lig 3/4 
; 7 - 12 0/7 S77 
i 13. = 19 0/7 2/7 
0 4/7 4/7 
- 1 - 6 3/7 D7 7 
‘ 7 - 12 0/7 Ly] 
: 13 - 19 1/7 a 


D. RECOMMENDATIONS FOR FURTHER WORK 
To reduce analysis computing time, the HP-9825 should be 
equipped with additional memory capacity. To eliminate use 
of tape for data storage, memory should be increased to 25K. 
The data analysis for each run should be stored and 
analyzed for mean and standard deviation. Several runs on 


the same data set should be performed to verify distribution. 


a7 


Number of scans sampled in each run should be increased 
by either solving signal trigger problem at low tape speed 
or by reprogramming NIC-80 computer to perform and store 


partial analysis results. 
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APPENDIX A 


Operating Instructions for Data Analysis 


A. Initialization 
1. NIC-80, Analysis program not loaded 

(a) Turn on computer with key, depress ''STOP" 

(b) Turn on Data Bus unit 

(c) Set toggle switches 0 - 1l to up position 

(d) Depress "LOAD PC" 

(e) Depress "CONTINUE" 

(f) Press "EXECUTE" 

2. HP-9825 

(a) Turn on HP-9825 

(b) Insert Data Analysis cassette and load track 0, 
file 1 for NIC-80 loading program 

Ce) Push “RUN” + -¢; dy 15 Fy 1) Oy Ly S, © is displayed. 

(d) Push "n'" and "CONTINUE"; load file # is displayed. 
The NIC-80 data analysis program is located on files 2, 3, 4 
and 5. Enter one of these numbers and push ''CONTINUE". 

(e) Display =e, d, 2,2, ny, 0; 1, Ss, t. Push "1" 
and continue. The section of the program will be transferred 
to NIC-80. At the end of each transfer, check sums are 
printed and should be identical. 

(£) Repeat steps (d) and (e) until entire program is 
transferred. If NIC-80 contained MIF analysis program pre- 


viously, only files 2 and 3 are necessary to transfer. 
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(g) On NIC-80, push "STOP", depress ''START" and push 
"EXECUTE" on completion of program loading. 
3. Tape Recorder 
(a) Turn switch to 3.75 IPS 
(b) Switch on DC bias boxes. These should be turned 
off during non-operating periods to avoid draining batteries. 
(c) Set frequency control knob on appropriate channels 
to FM and lke bandpass for signal and trigger. Voice channel 
requires 37.5 IPS and 100 Ke for operation. 
(d) Turn on oscilloscope. 
4. Pulse Generator 
(a) Turn on power 
(b) Set to external trigger pulse. 
(c) Set width to 10 microseconds. On tape recorder, 


depress "DRIVE" and adjust oscilloscope for optimum display. 


B. EXECUTION 
1. HP-9825 
(a) Load track 0, file 8 and push "RUN" 
(b) Data entry; 
Wavelength, e.g., .6328 
Optics, e.g., .4572 
Dwelltime, e.g., 95 
Number of scans, e.g., 125 
| Push "CONTINUE" after each entry. On last ''CONTINUE' NIC-80 
| will enter Display/Standby mode. HP-9825 will display 
"CONTINUE WHEN READY". 
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2. NIC-80 

(a) Verify green control light is on 

(b) Adjust CRT for optimum display of scans 

(c) Adjust DC level by knobs on bias boxes. Use low 
resolution on vertical display scale, and adjust signal to 
center around horizontal scale. 

(d) Adjust horizontal position of signal by pulse 
generator delay knob so as to center around vertical axis. 

3. To start data collection, push ''CONTINUE" on HP-9825. 
When all scans are completed, NIC-80 will return to Display/ 
Standby mode and may be taken off-line. 

4. Data Analysis (HP-9825 only) 

(a) On completion of data collection, "HISTOGRAM AND 
CONTINUE" is displayed. Check that plotter and printer are 
ready before pushing ''CONTINUE". 

(b) Data entry; 

Date, e.g., 021577 (month, day, year) 
Time, e.g., 1423 
Cn arg, e.g., 8.07 
Power, e.g., -8 
Push "CONTINUE" after each entry. 

(c) Unnormalized values for MIF spatial frequencies 0 
and 1 are printed for verification of absolute and relative 
values. 

(d) Data entry; "No. of class intervals" is displayed. 


Enter number of histogram bins per unit standard deviation. 
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Normally 3. Push ''CONTINUE". 
(e) Display ''FREQ PLOT, YES-1". Push "1" and '"'CON- 
TINUE" if plot of MIF mean + 1 standard deivation for spatial 
frequency 0 - 19 is required. 
(f) Display "LN DATA? YES-1". Enter 1 and push 
"CONTINUE" if data analysis of normal distribution is completed. 
Display "LN ANALYSIS". Push "CONTINUE" if data 
is to be converted to Ln data. 
(g) Display "DATA SUM? YES-1"'. Enter 1 and push 
"CONTINUE" if data summary is required. 
Display ''SUMMARY, CONT.". Push "CONTINUE" to 
start execution of data summary program. 
(h) Display "SPAT FREQ". Enter spatial frequency 
number (0 - 19) for which histogram will be plotted. 
Display ''REPLACE PAPER AND CONTINUE". Push 


"CONTINUE" when plotter is ready. 
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if X=l;gto "START" 

mdec;par l;fmt 1,fz10.5;enp "wavel",W,"opt",0O 

enp “DWELLTIME",r0,"# OF SCANS",N 

20>+P;dim F[P],A({P],M[P],S[(P],X[(P];wrt 2,1,W,9,r0,N;par 0 
for J=l to 20;rdb(2)+A;rdb(2)+B;512*B+A>+F [J] 

if F[(J]>524287;F[(J]-108576+F([J] 

next J;trk 1 

enp “CONTINUE WHEN READY",r0 

wtb 2,r0 

for I=l to N;for J=l to 20;rdb(2)+A;rdb(2)+B;512*B+A+A [J] 
10: if A[{J]>5242387;A[(J]-1048576+A[J] 

ll: rdb(2)+A;rdb(2)+B;512*3B+A+r0 

12: 1£ r0>5242387;r0-1048576-+r0 

13: A[J]*A[J]+r0*r0+r0;/r0+A[J] 

14: A[J])/F[J]+A[J] ;A[J]*A[J]+S[J]+S[(J] sM[J]+A[T]>M[T] 
15: next J;rcf I,A[*] 

16: next I 

17: dsp “HISTOGRAM ?,and continue";stp ;flt 5 

138: enp "DATE",r99,"TIME",r98,"CN ARG",r94,"POWER",U 

19: "PLOT" :for Q=l to 20;trk 1 

20: (S[Q]-M{Q]*M[Q]/N)/(N-1)+S[O] -M[Q] AN*M[O]-¥S[O] +S [9] 
21: next Q 

22: flt S:prt "MTFO" ,M[1],"MTF1L",M[2];if r97=l;gto "DS" 
23: enp "NO CLASS INTERV",T;beep 

24: "START" :enp "FREQ PLOT, YES=1",r92;if r92=l;gto "“FREPLO" 
25: if r97=l;gto "DS" 

26: "L"senp "LN DATA?, YES=1",r97;if r97=1l;l+Z;gto "LN" 
27: "DS":0+X;enp "DATA SUM?,YES=1",X;if xX=l;trk O;ldf 9 
28: "PLOTO":enp "SPAT FREQ",A;A+1+Q;trk 1 

29: asp “REPLACE PLOT PAPER AND CONTINUE";beep;stp 

30's MlOl=3*S (OleCyM{Ol4+3*5S (CO) +. 

31: S[Q]/T+R;C-2,2*R+r9l;scsiz 1,6,1,5,1 

32: scl C-2*R,D,-.2,1.7;plt D,0,l;plt C-R/10,0,2;fxd 1 
33: plt r91+R/2,0,l;plt r91+R/2,,7,2;for I=0 to .7 by ol 
34: ltr r91+R/1,5,1I;lbl I;plt r91+R/2-R8/10,1,1 

35: plt r91+R/2+R/10,1I,2;next I 

36: -3*T+r88;fxd 0 

37: for I=C to D by R;olt I,.01,l;plt I,-,01,2 

38: ltr I-R/3,-.08;lbl r88 

39: r68+l>r388;next I;pen;csiz 1.6,1,5,1 

40: ltr C+3*R,1,6;lbl "RELATIVE DISTRIBUTION OF MTF" 
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Figure 8A. Data Analysis Program Listing 
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pen;plt M[Q],0;plt M[Q],.7,2 

ltr C+R,-.15;lbl “UNITS = 1/",T," STANDARD DEVIATIONS" 
"CURVE":for I=C to D by (D-C)/100;(I-M[Q])/S[Q]+rl 
cl*rl1/2*rl;exp(-rl)*R/(S[Q]*y¥ (2*7) )+r2 

plt I,r2,-2;next I;C+r2;r2-R-r3 


: Cer2s;r2+Rer3;f1lt 5;A+A 


for K=1l to 11+T*6;0+rK;next K 

"HIST":for I=l to N;C+r2:r2+R-r3 

ldf I,A[*];for K=l1l to 11+T*6 

flt 5S;1if A[Q]>r2;if A[Q]<r3;rkK+l>rkK 
r3*r2;r3+R>r3;next K;next I 

Cer2;r2+Rer3;for K=1ll to L1+T*6;rKerl 

rl/Nerl;pen;plt r2,0,-2;0lt r2,rl,2;plt r3,r1,2 

plt r3,0,2 

(r2-M[0))/S(Q]*r7; 27*r7/2>r7 2exp (-27) *R/ (S101 *7 (240) ) +44 
(rc 3-M[Q])/S[Q]*r7; r7*r7/2er7 3 exp (-r7) *R/(S[Q]* 7 (2*7))er5 
R*r4t+abs (r5-r4)/2*R>r8;r8/R>r8 

(N*¥*r8-rl*N)* 2/(N*r8)er9;r9+r1l0>rld 

r3>r2;r3+Rer3;next K 

"LABEL" :C+R/3erl;fxd O;ltr rl,1,.45 

lbl "SPATIAL FREQUENCY",Q-1 

1£ cr97=l;ltr r1l,1,25;lbl "TEST DIST LN NORMAL";gto +2 
ltr r1l,1,25;lbl "TEST DIST NORMAL" 

ltr rl,1.35;1bl "NO OF SCANS",N 

fxd 4;ltr rl,1.15;lbl "WAVELENGTH",W;: ltr r1,1,05 

lbl "OPTICS",O 

fxd O;ltr rl,.9;1lbl "DATE",r99," TIME", r98;flt 5 


: M{OQ] +R/3erl; itr ri,i»35sibl "St, dev",Sto] 


ltr r1l,1,45;lbl “MEAN",M[Q] 

Ltr rvrl1,1.25;lbl "CHI SQ",rlO;ltr cl,1.15;f£xd 3 
lbl "CN",r94;f£xd 0 

IbL. *E";U 

for I=l to 90;0+rI;next I 

if r97=l;gto "Ds" 


geo 8h" 


"LN":dsp "LN ANALYSIS";stp ;trk Ll 

tor J=l to 20;0+M[J]+S[J];next J 

for I=l to N 

ldf I,A(*];for J=l to 20;if A[J]=0;if J=20;gto +5 
1£ A[J]=O;next J 


Figure 8B. Data Analysis Program Listing 
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In(A{J])+A[J] 


: A(J]*A(J]+S[J]+S[J] sM[J]+A[J] +M[J] 


next J;rcf I,A[*] 

next I;gto "PLOT" 

"FREPLO":scl -3,20,-,25,1.5;fxd O;csiz 1,6,1.5,1 

plt 20,0,l;plt 0,0,2;plt -1.385,0,l;plt -1.85,1,2 

for I=l to 20;plt I-1,.01,l;p1lt I-1,-,01,2;next I 

for I=l to 20 by 5;I-l>r95;ltr r95-,3,-.,05;lbl 195 
next I;fxd l 

for I=0 to l by ,l;ltr -1,85,I;lbl I;plt -1,9,I1,1 

for 1=0 to 1 by .lsltr -1,85,13;1bl I;plt =—159,1,1 

olt -1,32,1I,2;next I 

ltr 2,-,.19;lbl “SPATIAL FREQUENCY IN 1/ MILLIRADIANS" 
ltr 2,1,45;lbl "MTF MEAN AND + - 1 STANDARD DEVIATION" 
fxd O;sltr 2,1,35 

lbl "DATE",r99," TIME",r98," NO OF SCANS",N 

fxd 4:ltr 5,1,25;lbl “WAVELENGTH",W," OPTICS",0O 

flt Ssltr 5,1,.15;fxd 3;lbl "CN",r94;fxd O;lbl "E",U 


: M[lL]+r93 


csiz 3,1,1,90;ltr -2.1,,4;lbl "MTF" 

for Q=l1 to 20;50+Q4+r48;70+0+r493M[0] /r93+r4d 

S[Q] /r93+r49 

plt Q-1,r438-r49,l;plt G@1,r48+r49,2 

plt Q-1.05,r4d,l;plt Q-.95,r48,2;pen;next Q;gto "L" 


Figure 8C. Data Analysis Program Listing 
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dsp "SUMMARY ,CONT";beep;stp 
"X"sfor Q=lL to 20 
$[0]/T+R 
id [QO] -3*510)]>C;M(Q]+3*S(Q]+D 
for K=1ll to 11+T*6;0%rK;next K 
"HIST":for I=l to N; Cor2;r2+Reor3 
trk Ll 
ldf I,A[(*];for K=lLl to 11+T*6 
flt 5;if A[Q]>r2;if A[Q]<r3;rkK+lerk 
r3er2;r3+Rer3;next K;next I 
Cer2;r2+R>r3;for K=1ll to Ll+T*o;rkerl 
cl/Nerl 
(r2-M[Q ea r7*r7/2er7;exo(-r7) *R/(S[Q] 
(c3-M[OJ])/S[OJ]er7;3 r7*r7/2er7 3 exo(-r7) *R/(S[Q]} 
R*r4+aps(r5-r4)/2*R>r38; r8/Rer8 
(N¥r8—r1l*N)* 2/(N*r8)%r9;r9+r1l0+rld 


)>r4 
{ei 


r3er2;r3+Rer3;next K;rl0+X[Q] ;0+rl0;next Q;gto "PRINT" 
"PRINT" :fmt Z,25x,"DATE",2x,£6,0,2x,"TINE",2x,£4.0 


fmt 1, 25%, "Cn", 2% i523 ek pS" £3 60 p 2h" OPTICS’ 7 


fmt 3,35x,"DATA SUMMARY" 
if Z=l;gto "Ll" 
wrt 6,3,13 

wtb 6,10,10 

wrt 6.2,r99,r98 
wtb 6,10,10 

wrt o,1,r94,U,0 
wtb 6,10,10 


fmt 4,25x,"TEST DIST.NORMAL",26x,"TEST DIST LN NORMAL" 

wrt 0.4 

wtb 6,10,10 

fmt 5,15x,"SPAT FR",5x,"MEAN",1lL1x,"ST,DEV",10x,"CHI SQ",10x, 


wrt 6,5;wtbd 6,10,10 
fmt 615%, 206%, 60. oles 1 tce Sect. 


for I=l to 20;wrt 6,.6,I-1,M[1I],S[I],X[I];next I 


for I=l to 2U0;wtb 6,27,10;next I 
fmt 7,72x,ell,5 

if Z=l;gto "Ll" 

trk O;ldaf 3 

"Ll"sfmt 7,72x,ell,5 

for I=l to 20;wrt 6,7,X[I];next I 
trk U;ldf 8 


Figure 9. Data Summary Program Listing 
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2x,£6.,4,2x 


o« 


: ra 
a ogi 


110600 
4304 
44034 
10377 
402340 
0 
2175677 
4304 
3001737 
3001745 
16016 
3001737 
3001745 
1614 
3001737 
3001745 
1600 
3001737 
3001745 


2703677 
27 
3001742 
1600 
3001715 


Figure LOA. 


MEMA(600 /READOUT LIGHTS 
LCWORD /CONTROL WORD 
STU /READ STATUS 


ANDA(377 /COMPARE SETTINGS 
A-MZ(340 /CHECK HARD SETTING 
JMP 0 /RETURN CHECK LOOP 


ZERMA 1677 
LCW (0 
JMS@INP 
JMS@1745 /PUTAC 

1616 /ADR OF WAVELENGTH 
JMS@INP /INPUT 

JMS@1745 /PUTAC 

1614 /ADR OF OPTICS 

JMS@INP /INPUT 

JMS@1745 /PUTAC 

1600 /ADR OF DWELLTIME 
JMS@INP /INPUT 

JMS@1745 /PUTAC 

1620 /ADR NO CF SCANS 
MEMA(20 /COUNT VALUE 

ACCM 1677 /COUNT INDEX 

JMS 1500 /CAL ADRESS 
MEMA@1576 /CAL VALUE 
ANDA(777 /9 BITS 

HSP /IS HP 9825 READY 

JMP 31 /WAIT LOOP 

PHS /WRITE TO HP 9825 
MEMA@1576 /CAL VAL RECALLED 
RIS ll /RIGHT SHIFT 

HSP /IS HP 9825 READY 

JMP 36 /WAIT LOOP 

PHS /WRITE REMAINING DATA 
MPOM 1576 /INCR, CAL ADRESS 
MMOZ 1077 /DECR.&TEST 

JMP 27 /RETURN LOOP 


70 TO ACC & 1677 
/LOAD CONTROL WORD 
/INPUT 


JMS@1742 /GETAC 
1600 /ADR, OF DWELLTIME 
JMS@1715 /FIX 


Nic 80 Program Listing 
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’ 
tk 


™ 


a > ey Ory 


3111722 


3001723 
2111704 
2405677 
2111626 
2405672 
2111712 
4304 
4302 
5220 
100 

0 

) 
110200 
4304 
2111700 
2405677 
2711701 
2405076 
110040 
2405675 
2165674 
3111677 
2505674 
3111676 
2505674 


MEMA@1722 /FACM 

LDWELL /LOAD DWELL TIME 
MEMA(OQ /0 TO ACC 

LDELAY /LOAD DELAY TIME 

MEMA 1711L /STAND BY LOOP 
LCWORD /LOAD CONTROL WORD 
SET /SET MEASURE FF 

STOP /READ 1 SCAN 

HSRF /SKIP ON HP REQUEST 

JMP 53 /RETURN STAND BY LOOP 
RHSR /READ HP 

A-MZ(314 /L FOR FUTURE BRANCH 
JMP 64 /FUTURE BR.ADR, 
JMS@1723 /SHORT TERM RESET 
MEMA 1704 /B4(0) 

ACCM 1677 / 

MEMA 1626 /# OF SCANS 

ACCM 1672 / 

MEMA 1712 /CONTROL WORD #2 
LCWORD /LOAD CONT,WORD 

SETM /SET MEASURE FF 

STOP /READ 1 SCAN 

JMP 100 / 

0 /NOT USED 

QO /NOT USED 

MEMA(200 / 

LCWORD /DISPLAY CONTROL 

MEMA 1700 /8B0O(0),LEFT END 
ACCM 1677 / 

MMOA 170L /80(1777),RIGHT END 
ACCM 1676 / 

MEMA(40 /32 POINTS 

ACCM 16075 / 

ZERM 1674 / 

MEMA@16077 /INTEGRAL BASE LINE 


A+MM 1674 / 
MEMA@1676 = = / 
A+MM 1074 / 


Figure 10B. Nic 80 Program Listing 
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2125677 
2705676 


2463701 
142 
3001776 


2703626 
71 
53 


Figure LOC. 


MPOM 1677 / 
MMOM 1676 / 
MMOMZ 1675 / 
JMP lll / 

MEMA 1674 / 


RASH6 /DIVIDE BY 6 

ACCM 1674 /AVERAGE DC COMPONENT 
MEMA 1700 /BO(0) 

ACCM 1677 / 

MEMA@1677 / 

A-MA 1674 /SUBTRACTING 
ACCM@1677 


/ 
MPOMA 1677 / 
A-MZ 1701 / 


JMP 126 /RETURN ADR,LOOP 
MEMA( 2 

JMS@1722 /MOVE BO+B2 
MEMA 1703 /B3 ZERO LOOP 
ACCM 1677 / 

MEMA 1700 /80 


ACCM 1676 / 
ZERM@1676 / 
ZERM@1677 / 
MPOM 1677 / 
MPOMA 1676 / 
A-MZ 1701 / 


JMP 142 /RETURN ADR, LOOP 
JMS@1776 / 
JMS@1622 /DFPT DATA TRANSFER 


ASCNTR /ADVANCE SWEEP COUNTER 
MMOZ 1626 /DECREASE SCAN COUNT 
JMP 71 /RETURN ADR,.NEXT SCAN 

JMP 53 /RETURN ADR.STAND BY L,. 
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152 
110024 
2405675 
2111702 
2405677 
2111703 
2405676 
3111677 
10777 
0474 
1456 
4473 
3111677 
405031 
6474 
1463 
4473 
3111676 
10777 
6474 
1470 
4473 
3111676 
405031 
6474 
1475 
4473 
2125677 
2125676 
2707675 
1454 
1001445 


Figure l1OD. 


152 /RETURN ADRESS 


MEMA(24 /# OF SPAT.FR, 
ACCM 1675 /CONTING INDEX 
MEMA 1702 /B2,REALS 

ACCM 1677 /CURRENT 32 
MEMA 1703 /B3,IMAGINARIES 
ACCM 1676 /CURRENT B3 
MEMA@1677 /VALUE B2 
ANDA(777 /9 BITS 

HSP /HP READY 

JAP 1456 /WAIT LOOP 

PHS /WRITE TO HP 
MEMA@1677 /RECALL 82 

RIS ll /RIGHT SHIFT 

HSP / 

JMP 1463 / 

PHS / 


MEMA@1676 / 
ANDA(777 / 


HSP / 

JMP 1470 / 

PHS / 

MEMA@1676 / 

RIS 1l / 

HSP / 

JMP 1475 / 

PHS / 

MPOM 1677 /INCREASE 82 
MPOM 1676 /INCREASE B3 
MMOMZ 1675 /DECREASE & TEST 
1454 /RETURN ADR, 


JMP@1445 /RETURN TO MAIN PR, 
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